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Abstract. The metal-insulator transition and the problem of d-electron delocalization are investigated 
in the pyrite system NiS2-a;Sea; under pressure using density functional theory (DFT). We test several 
approximations, including the generalized gradient approximation (GGA), the GGA+U approach, and 
hybrid functionals. In addition we apply the GW approximation and perform Hartree-Fock calculations. 
The important role of the chalcogen dimers in the electronic structure is discussed within GGA, which 
sufficiently describes the role of the lattice in the metal-insulator transition. In addition, the magnetic 
phase diagram is determined. However, the electronic properties are inadequately described - the insulating 
ground state of NiS2 cannot be obtained - and the magnetic order is slightly overestimated. If correlations 
are taken into account within GGA+U, the insulator is found, but the non-magnetic ground state of 
the doped samples is not accessible. Mixing Fock exchange with local approximations as GGA correctly 
reproduce the insulating ground state in NiS2, and metallic ground states with doping and pressure. The 
insulator is - in contrast to earlier suggestions - of Mott type, i. e. the gap opens between Ni 3deg states. 
Due to the Fock term, the magnetic order is strongly overestimated. The applied GW calculations are not 
able to correct the metallic character of NiS2 and underestimate the gap by one order of magnitude. 

PACS. 71. 15. Mb Density functional theory, local density approximation, gradient and other corrections 
- 71. 20. Be Electron density of states and band structure of transition metals - 71.30.-|-h Metal-insulator 
transitions and other electronic transitions 



2 Cosima Schuster: Electronic and magnetic properties of NiS2-i:Sei:: a comparative study 

1 Introduction electronic structure is explained by filling the band struc- 



Crystal field splitting, electron-electron interaction, and 
the 3d band width are almost equal in the pyrites MX2 
(M=Fe,Co,Ni,Cu,Zn; X=S,Se) which gives rise to a large 
variety of electrical, magnetic and optical properties in 
these compounds. Of special interest is the metal-insulator 
transition in NiS2 with Se-doping [1] or under pressure 
[2]. The transition is not accompanied by a change in the 
lattice symmetry and is commonly believed to be driven 
by the electron-electron interaction. In addition, a low- 
temperature antiferromagnetic phase is observed both in 
the insulating and metallic regime of NiS2_a;Se2; [Tj and 
by applying pressure in NiS2 [3 . The electronic and op- 
tical properties of NiS2 were first investigated already 30 
years ago and re-investigated experimentally in the 
late nineties |5I6I7I9I8] . The renewed interest in transition- 
metal oxides and their metal-insulator transitions traces 
back to the attempt to understand the phase diagram of 
the high-Tc superconducting perovskites. Moreover, the 
attempt to use the metal-insulator transition of V2O3 or 
VO2 for technical applications triggered many studies up Hence, we clarify the role of the S-S dimers and dis- 



ture with the rf-clcctrons of the transition metal by going 
from Fe {3d^4s^ ) to Zn {3d^^4a'^). The band insulators are 
described quite well within the local density approxima- 
tion, however the microscopic origin of the insulating na- 
ture of NiS2 is not well understood. Since more advanced 
approximations within DPT are available now, several at- 
tempts to improve the material specific theory have been 
performed. Recently, Perucchi et at |13| have discussed 
two different mechanism for the metal-insulator transi- 
tions based on LDA calculations. Compression and expan- 
sion of the lattice alter the metallic state, albeit through 
two different microscopic mechanisms. On the other hand, 
Kunes et al. [H], using LDA+DMPT, have traced back 
the metal-insulator transition to a single control parame- 
ter, namely the size of the p gap. They confirmed, how- 
ever, that compression and doping influences this control 
parameter in different ways. Both studies concentrate on 
the electronic structure in the paramagnetic phase and 
neglect the magnetism. 



to now [TU]. cuss the energy scales mentioned in [13114] . Por the cal- 

Theoretical attempts using density functional theory culations we use the wien2k code [15 and the general- 

and the local density approximation (LDA) to describe ized gradient approximation (GGA) in the parameteriza- 

the electronic structure of the pyrites date back to 1987 tion of PBE [TH]. In a next step we determine the mag- 

|11| . with the focus on PeS2, where the important role of netic phase diagram, Sec. III. Surprisingly, GGA describes 

the S(e)-S(e) dimers present in the pyrite structure has the magnetic phase diagram quite well, slightly overesti- 

been pointed out. The role of electron-lattice interaction mating the phase boundaries as compared to experimen- 

when applying pressure was studied ten years later, again tal data. In contrast to GGA, LDA always favors the 

in PeS2 yL2j. Other pyrites are mentioned in ll , and their non- magnetic ground state. Calculations using improved 
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exchange-correlations functionals are presented in Sec. IV. estimated in experiments 24 Ij. The first configuration 

We present results of GGA+J7(SIC) [17] and EECE hy- contains S2 and Se2 dimers, the second only S-Se pairs, 

brid [18ll9j calculations. To compare with we have per- Thereby the second configuration is favored by 0.015 Ryd 

formed also GW calculations on top of LDA using ABINIT when the volume is fixed. In general, S-S, S-Se, and Se-Se 

[20121] . We focus on the advantages and shortcomings of pairs are present in doped samples [8]. 
the different techniques with regard to correlated materi- 
als. 

Using the mentioned lattice constants the internal po- 
sitions are obtained by a force minimization. The distances 

2 The role of the dimers for the electronic are summarized in Tab.dJ In NiS2, the S2 distance is found 

. . to be about 2.1 A, slightly (2 %) larger than the with ex- 

structure ^ b j \ > & 

perimental data d— 2.07 [1118] . Applying moderate pres- 

The pyrite structure is a simple cubic lattice with space sure, the dimer decreases only slightly from 2.104 to 2.095, 

group Pa3. It is best described in terms of the NaCl struc- thus the dimer-distance is changed by 0.4% when chang- 

ture with the transition metal in one sub-lattice and the ing the lattice constant by 3.5%. Applying even more pres- 

center of mass of the chalcogenide pairs in the other. Thus, sure, corresponding to a 20% volume effect, it decreases 

the S/Se atoms form distorted octahedra around the Ni only to 2.07 A. Also when expanding the volume drasti- 

atoms. In the following, we use the experimental lattice cally, the dimer-distance vary hardly, compare [14] . For 

constant of NiS2, a — 5.687 A [22 . Applying hydrostatic this reason, the metal-insulator transition applying pres- 

pressure is simulated by a shortening of the lattice con- sure is related to a volume effect, since the dimer-distance 

stant - 5 %/10 % volume reduction correspond to 1.7 is nearly less constant when varying the volume. Turn- 

%/3.5 % reduction of the lattice constant and to 5/12 ing to the doped samples, we find a large variation of the 

GPa, respectively. In the attempt to understand the phase dimer-distances with doping. In NiSe2, the Se-Se distance 

diagram with doping we concentrate on the end points is given by d = 2.49 A, which is about 4% larger than 

with x = and a; = 2. In addition we consider also x = 1, the experimental [llj . Worth mentioning is that the Se- 

i. e. NiSSe. For NiSe2 we also use the experimental lattice Se distance is larger than the Ni-Se distance. Turning to 

constant, a = 5.960 A ;22J. In case of NiSSe, we study two NiSSe, we observe that in NiSSe I the the S-S/Se-Se dis- 

different configurations within the pyrite unit cell - recent tances are comparable to dimer-distances of the pristine 

STM measurement on doped samples show no sign of a materials, see Tab. reftab. In NiSSe II, on the other hand, 

superstructure [53] - and use the theoretically obtained the S-Se distance corresponds to the averaged value of 

volume, athco — 5.87 A which is about 1 % larger than (i=2.3A. 
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Table 1. Lattice constants and typical bond lengths - given 
in A- of the systems under consideration, NiS2, compressed 
NiS2, expanded NiS2, NiSSe, and NiSe2. In addition, the energy > 



40 



differences of the bonding and anti-bonding splitting of the m 

O 



chalcogen (S, Se) s and p states are given. 



System 


a (A) 


d: Ni-S/Se 


d: S(e)-S(e) 


NiS2 


5.28 


2.21 


2.07 


NiS2 


5.49 


2.31 


2.10 


NiS2 


5.69 (exp.) 


2.40 


2.10 


NiS2 


6.04 


2.56 


2.11 


NiSSe I 


5.87 


2.48 (S) 


2.14 






2.44 (Se) 


2.42 


NiSSe II 


5.87 


2.41 (S) 


2.30 






2.50 (Se) 




NiSe2 


5.96 


2.48 


2.49 


System 


a (A) 


AEs (eV) 


AEp (cV) 


NiS2 


5.28 


2.2 


11.0 


NiS2 


5.49 


3.2 


9.9 


NiS2 


5.69 (exp.) 


3.4 


9.0 


NiS2 


6.04 


2.3 


7.2 


NiSSe I 


5.87 


3.35 (S) 


8.6 (S) 






2.3 (Se) 


8.0 (Se) 


NiSSe II 


5.87 


1.9 


7.7 


NiSe2 


5.96 


1.8 


7.4 



Next, we study the volume/dimer-distance effect on 
the electronic structure, see Fig. [l] where we show the 
density of states (DOS) of NiS2. Generally, the presence 
of S/Se2 dimers in the pyrites leads to a splitting of the 
bonding and anti-bonding S 3s states as well as S 3p states, 
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Fig. 1. (Color online) Density of states of NiS2. We show the 
S 3s, S 3p3, and Ni 3d states. The bonding and antibonding 
parts of the S 3p states are indicated. 

as indicated. A change of the volume is expected to modify 
the bandwidth. With increasing volume the bandwidth 
should decrease, and vice versa. However, the situation in 
NiS2 is more complicated as pointed out for FeS2 in jl2j . 
Following the detailed analysis in [T^] the peaks at —16 
and —12 eV are assigned to S 3s states. These states are 
observed in the XPS valence spectra at —15 and —11 eV 
[Sj. Thus, the splitting is correctly described in GGA, but 
the position shifted by 1 eV. Dominant S 3p contributions 
are found at the lower band edge (at —7 eV) and the 
upper band edge (at 2 eV). They form in itself a double 
peak structure, compare H^. The Ni 3d eg states - located 
between —4.5 eV and —2.54 eV and between —1 eV to 1 
eV - overlap with S 3p states. The Ni 3d t2g states, located 
at —2 eV, show almost no overlap with the S 3p-states. 

Applying pressure, see Fig. [2l upper panel, where we 
show a comparison of the DOS of compressed and un- 
compressed NiS2, a rigid shift of the t2g by —0.5 eV due 
to the enhanced crystal field splitting is observed in the 



o 

Q 



60 



40 



20 



NiSa, a = 5.69 
NiSs, a = 5.49 




-16 -14 -12 -10 -8 -6-4-2024 
E - Ef (cV) 



Cosima Schuster: Electronic and magnetic properties of NiS2-2;Se2,: a comparative study 5 

Upon doping, see Fig. [21 lower panel, where we show a 
a comparison of the DOS of NiS2 and NiSe2, the bonding- 
antibonding splitting is reduced due to the increased Se-Se 
distance, both for the s and p states, see Tab. [TJ as com- 
pared to NiS2. In NiSSe I, the low energy double peak 
structure due to the s states splits into four peaks due to 
the presence of both S-S and Se-Se dimers with different 
dimer lengths. This could be an experimental signature 
whether a structure with both S-S and Se-Se dimers is 
realized or only S-Se pairs are present. The dip at Ep is 
more pronounced in the NiSSe II than in the NiSSe I struc- 
ture. In addition, the conduction band width is smaller in 
structure II (taken at the same volume). To summarize, 
the dominant effect of doping on the DOS/band width is 
the decrease of the splitting of the S 3p-states and the re- 
sulting stronger overlap of the Ni 3d with the S p states, as 
commonly suggested. The stronger hybridization is identi- 
fied as the driving mechanism of the metal-insulator tran- 




-16 -14 -12 -10 -8 -6-4-2024 
E - Ef (cV) 



Fig. 2. (Color online) Upper panel: Density of states of NiS2 
at the experimental lattice constant and compressed NiS2. The 
straight (red) line corresponds to zero pressure, the dotted 
(green) line to 12 GPa. Lower panel: Density of states of NiS2 ^i^i*-*^ upon doping, 
(red line) in comparison to the density of states of NiSe2 (green 
line) . 



Up to now, we have analyzed the density of states 
over the whole band width. Turning to the underlying 
band structure, we take closer look at the states near the 



DOS. Due to this shift the Ni t2g- and the S p-states over- Fermi level. Interestingly, the changes in the dispersion 
lap. The states below —3 eV are subject to an additional £{k) around Ep are similar when applying pressure and 
shift, so that the lower band edge is about 0.6 eV lower with doping. For comparison, we plot the according band 
for a = 5.49 Athan for a — 5.69 A. Hence the overall band structures in Fig. [31 The changes in the band structure 
width is broader and the bonding anti-bonding splitting near the Fermi level with pressure can be focused on two 
increases when compressing the lattice. In addition to the symmetry points: The band at F (0.3 eV) moves down- 
shift, the binding S 3p-states show a slight broadening, wards, almost touching or crossing the Fermi level, respec- 



but the rigid band shift is the dominant effect. 



tively. Second, the four degenerate bands at M are shifted 
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Fig. 3. (Color online) Band structure near Ep of - from top 
to bottom - NiS2, NiS2 with 10% volume reduction and for 
NiSez. 



towards the Fermi level. In NiSe2, the band at R moves 
from the Fermi level to about 0.4 eV above. The hybridiza- 
tion due to band shifts or decreased bonding-antibonding 
splitting does not manifest itself near Ep. 



properties of NiS2-i:Sei:: a comparative study 

3 Magnetic phase diagram 

The phase boundaries of the AFM ordering of type I 
with pressure and doping can be reproduced within the 
DFT calculations using GGA. From the experience with 
LDA/GGA studies of other magnetic materials it is ex- 
pected that LDA underestimates and GGA overestimates 
magnetic order and magnetic moments, see for example 
[25j . In NiS2, we consider an antiferromagnetic alignment 
of type A, i. e. the magnetic structure consists of fer- 
romagnetic planes which are coupled antiferromagneti- 
cally. Using this structure no supercell is needed, in ac- 
cordance with 3 . The magnetic properties - calculated 
within GGA - are summarized in Tab. [21 The given mag- 
netic structure is stable within the GGA calculations, frus- 
tration effects play a minor role. As expected, NiS2 turns 
out to be a high-spin compound with /i = 0.73/iB. With 
increasing pressure, the magnetic moment as well as the 
energy difference to the non-magnetic case decreases. In- 
creasing the pressure up to 12 GPa, the antiferromagnetic 
ground state is still favored, but the energy difference to 
the non-magnetic solution tends to zero. NiSe2 turns out 
to be non-magnetic, with zero moment. At the critical 
doping (NiSSe), we find also an antiferromagnetic order, 
which can be suppressed by applying pressure. Unfortu- 
nately, a lower boundary of the magnetic transition by 
using LDA could not be determined since LDA favors the 
non- magnetic solution in all cases. This can be related to 
a volume effect. Within LDA the ground state volumes 
are much smaller than the experimental or GGA volumes 
which may have an impact on the spin-polarized calcula- 
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tions. For completeness, the data for the energy differences 
using LDA are also given in Tab. [2j 

The electronic properties change due to the magnetic 
order. Band structure and density of states are heavily in- 
fluenced by the magnetism, see Fig. 31 where we show DOS 
and band structure of antiferromagnetic NiS2. The degen- 
eracy of the bands is lifted from a fourfold to a twofold 
degeneracy. States near Ep, especially at R and M, are 
emptied. As a result, the dip in the DOS at the Fermi 
level becomes more pronounced. 



NiS2: antiferromagnetic spin alignement 




> 



o 
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Fig. 4. (Color online) Band structure and density of states of 
NiSa with anti-ferromagnetic spin configuration. 



Table 2. Magnetic ground state and magnetic moment (in ^b) 
within LDA/GGA calculations for the systems under investi- 
gation. AE — (-Enon-poi — -Bmag)/ (number of Ni) is the energy 
gain due to the magnetic ordering per Ni atom and is given in 
meV. 



compound 


Vxc 


mag. moment 


AE 


NiSa, 5.69 A 


GGA 


±0.73 


2.1 


NiSa, 5.69 A 


LDA 


±0.20 


-0.14 


NiSa, 5.49 A 


GGA 


±0.45 


0.17 


NiSa, 5.49 A 


LDA 


±0.11 


-0.17 


NiSa, 5.28 A 


GGA 







NiSSe, 5.69 


GGA 


±0.4 


0.034 


NiSSe, 5.87 


GGA 


±0.5 - ±0.6 


0.88 


NiSSe, 5.87 


LDA 


±0.15 


-0.14 


NiSea, 5.96 


GGA 








4 GGA+U, Fock exchange 

Local approximations to DFT like LDA and GGA fail to 
describe the insulator [13j. Recently, a variety of insulat- 
ing transition-metal oxides have been investigated using 
LDA+f/ 17J, LDA±DMFT [Mj, Fock exchange or hy- 
brid functionals [19] . or the GW approximation 26J. In 
this work, we apply GGA+C/, hybrid functionals, Hartree- 
Fock, and GW calculations to NiS2. We use the imple- 
mentation of GGA±J7, comparing the fully localized limit 
(SIC) as applied to NiO [17] and the around mean-field 
treatment (AFM) j27| for the double counting correction, 
and the EECE treatment [19] within the wien2k code [15] . 

To start with, we perform constraint GGA calculation 
to obtain a reasonable value for U. Following the treat- 
ment in |28j and assuming a Ni d^ configuration as in 
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NiO, we obtain U^s = U - J = 6.39 eV for NiSa. C/cff 
decreases to 5.89 eV for NiSSe, and further to 4.94 eV 
for NiSe2. The different chemistry such leads to a con- 
siderable decrease of the correlation energy with doping. 
With decreasing volume, [/off is constant. In the following, 
we always use the antiferromagnetic spin structure as dis- 
cussed in the previous section. The DOS, calculated with 
UcS = 6.4 eV and the SIC correction, is plotted in Fig. [51 
upper panel, and shows a rather poor agreement with ex- 
perimental data [i6 7 . The Ni 3c? states are shifted to too 
low energies, forming a triple peak structure at the lower 
band edge. The S 3p states are homogeneously spread over 
the whole energy range. Using the experimental value of 
U = 5 eV the gap values, namely A^^'^ = 0.6 eV and 
^AFM _ Q 2 g^j.g jj-^ -j-j^g range of the experimental val- 
ues, but the Ni 3d states are located at —4 eV instead 
of —2 eV. A similar behavior was found in case of the 
ferromagnetic metallic pyrite C0S2 

While GGA+U is successful for insulators, the descrip- 
tion of correlated metals where the requirement of an in- 
teger orbital occupation and well established long range 
order [30] is no longer fulfilled may fail. In case of NiS2 un- 
der pressure, GGA+U correctly shows that this system is 
a non-magnetic metal. The results for NiSe2 on the other 
hand, are very sensitive to the choice of the method. Using 
SIC the ground state is a ferromagnetic metal, in case of 
AFM the antiferromagnetic metal is nearly degenerate to 
a ferromagnetic metal with two different moments. 

In a next step we perform calculation using Pock ex- 
change for the Ni within the muffin tin spheres (exact ex- 
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Fig. 5. (Color online) Total density of states of NiS2 (red line) 
and contributions of the Ni 3d (green and blue line) and S 
3p (pink line) orbitals: we compare the GGA-|-?7 calculation 
(left-hand side) with the hybrid functional with 20% EECE 
(right-hand side). 

change for correlated electrons, EECE) and LDA for the 
correlation potential. This approach is known to give quite 
reasonable results for correlated electron systems |18I19) . 
Instead of a parameter U and the choice for the double 
counting corrections which give deviating results for insu- 
lating NiS2 and metallic NiSe2, the mixing parameter of 
the hybrid functional has to be determined. As expected 
100 % EECE overestimates the gap by a factor of 5. It 



Cosima Schuster: Electronic and magnetic properties of NiS2-2;Se2,: a comparative study 9 

also shifts the Ni 3d states to the lower band edge, i. e. to band has dominant S character. The S character of this 

— 10 to —8 eV. In this way, 100 % EECE is comparable to band shows up only at F, compare FeS2 [12] ■ According 

U — 10.7 eV. The gap is induced within the S p states, and to the band structure the metal- insulator transition under 

not as expected between the Ni 3d (Mott-Hubbard type) pressure is triggered by the band with Ni character which 

or Ni 3d and S 3p states (charge transfer type). If 25 % provides the states at the Fermi level between F and X. 

EECE is added to LDA the Ni 3d(2g states are found at With doping the metallic state is attributed to the states 

—4 eV. The DOS shape corresponds to the DOS obtained between X and M. In this way, a metallic ground state is 

for U ~ 5.4 eV (suggested U value from [H]). Decreasing obtained for NiSe2 due to the increased dimer-distance. 

the amount of EECE further, the gap closes below 20%; On the other hand, reducing the volume also closes the 

therefore we perform the following calculation with a hy- gap at constant dimer-distance. The results for the band 

brid functional where 20% EECE is mixed to LDA. Then, structure are shown in Fig. [SI where we compare insulat- 

gap in the DOS opens within the Ni-S hybridized states, ing NiS2 with metallic NiS2 under pressure and metallic 

located between —2 and 1 eV, at the preformed dip in the NiSe2. 

GGA calculation, as can be seen in Fig. [51 The bonding- To summarize, GGA+C/ is not able to describe the sys- 

antibonding S 3p states dominate at the lower valence and tem properly, neither with the calculated nor estimated \J 

higher conduction band edge. The double peak structure values. Hybrid functionals using Fock exchange can re- 

of the empty states observed in XAS can be traced back produce the insulator and give a reasonable DOS. This 

to a Ni 3deg double peak and the S 3p antibonding double approach fails, however, in describing the magnetic phase 

peak. Using this hybrid functional, an antiferromagnetic diagram by overestimating the ordered phase. Thus, we 

solution is obtained, although LDA itself favors the non- have identified a system where the advantage of hybrid 

magnetic ground state. The calculated moments become functionals versus GGA+J7 can be clearly demonstrated, 
smaller under pressure and with doping but do not vanish. 

5 Quasi-particle picture: HF and GW 

The gap of NiS2 is very sensitive to the S-S dimer- 
distance d. Using a structure with d = 2.15 Athe band Corrections to the energy eigenvalues in DFT can be esti- 
gap amounts to 0.3 eV, whereas it is 0.6 eV with d = 2.10 mated by first order many-body perturbation theory with 
A. Increasing the dimer distance to d = 2.36 A, i.e. using respect to the difference in exchange correlation self-energy 
the structural parameter of NiSe2 but the lattice constant and exchange correlation potential of DFT. Describing the 
of NiS2, the gap closes. A closer look at the weighted dielectric properties of semiconductors, the GW approxi- 
band structure shows that the gap is located at the mation has become the method of choice [31 . Thereby, the 
lower band has dominant Ni character, whereas the upper non-local, energy-dependent self-energy is approximated 
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NiS2 0.2EECE 



S 0.0 

c 

ai 



NiS2, a = 5.49 A 



S\ 0.0 



NiSe2 



is described correctly [26J. Even properties of V2O3 can 
be described [33 . In NiS2, the electron-lattice coupling 
is strong and found to be the driving mechanism of the 
metal-insulator transition, thus it can be expected that the 
GW approach is able to provide a description of the insu- 
lating state. Similar to VO2, in the case of NiS2, a self con- 
sistent scheme has to be applied to get rid of the metallic 
starting point. In the following we use the GW implemen- 
tation of the ABINIT package [21]. We compare, see [34], 
Hartree-Fock calculations with self-consistent, but static 
GW schemes. The self-energy can be written as a sum of 
the exchange self-energy (static screened exchange), the 
correlation self-energy (Coulomb-hole term) and the dy- 
namical (energy-dependent) self-energy, and the approx- 
imations used within GW are classified along this parti- 
tioning. Here, we compare the screened exchange with the 
screened exchange plus Coulomb-hole. Another static GW 
approximation is the quasi-particle self-consistent GW [3S] . 
As we start from the Kohn-Sham orbitals obtained within 
Fig. 6. Band structure obtained using a hybrid functional with lqa, spin-polarized calculations are in accordance with 





20% EECE. Top: NiS2; middle: NiS2 under pressure; bottom: 
NiSe2. 

by the Green's function and the screened Coulomb in- 
teraction, T, = GW 32\ The screening is calculated on 

the level of the random phase approximation. Within the 

(3/8,3/8,1/8) between X and M. As it can be assumed 
"one-shot" procedure Gq and Wq are calculated from LDA/GGA, 

that the correction do not vary in the vicinity of the cho- 

but this Go Wo approach gives appropriate results only if 

sen k-points, Kl should give the correction at K2 at X, 

GW is a small correction to LDA. Transition-metal ox- 

K3 at M, and K4 at R. 

ides can be treated rather well using self-consistent ap- 



non-polarized calculations since LDA favors a non- magnetic 
ground state. For simplicity we chose the k-points where 
the corrections are calculated out of the optimal k-set. The 
given k-points Kl (1/8,1/8,1/8) and K4 (3/8,3/8,3/8) are 
located between F and R, and K2 (3/8,1/8,1/8) and K3 



proaches. For example, the electronic structure of VO2 
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Table 3. Energy gaps obtained within diflterent approxima- 
tion used for the self consistency, namely screened exchange, 
screened exchange plus Coulomb hole, and model GVF[35]. 



method 


spin 




k-point 


gap 


HP 


non- 


■pol 


r 


0.3 eV 


HF 


pol 




r 


0.4 eV 


scrEX 


pol 




M 


0.175 eV 


COHSEX 


non- 


■pol 


M 


0.021 eV 


COHSEX 


pol 




M 


0.025 eV 


QPscGW 


pol 




M 


0.027 eV 



On all GW levels, gaps are induced to the band struc- 
ture of NiS2 at the chosen k-points, where the smallest gap 
occurs at K3 (M). The results are summarized in Tab.[3l 

The behavior is in accordance with the GGA calcula- 
tions as we can assign the dip at Ep of the GGA DOS 
to the bands crossing at M. However, it is found that one 
band crosses the Fermi level between R and F. For this 
reason, GW calculations cannot describe the insulating 
state of NiS'2, even though a gap is formed at all given 
k-points. Therefore, calculations of compressed NiS2 and 
NiSe2 lead to very similar results for the gaps. Moreover, 
the metallic phases show larger gaps at these k-points, but 
the crossing is present in each case. 

On the other hand, the smallest gap is loacted at F 
within Hartree-Fock, in accordance with the results dis- 
cussed in the previous section. In addition, the gap values 
resulting from the calculations with 20 % EECE and from 
the Hartree-Fock calculation are astonishingly of similar 
magnitude. However, the total density of states shows no 



gap at Ep but a large gap between —5 to —1 eV, despite 
the direct gaps at the chosen /c-points. 

As a major shortcoming, the GW and Hartree-Fock 
calculations, on top of the non-magnetic starting point 
generated by the LDA, do not generate a magnetic ground 
state as the EECE does. The local magnetic minimum is 
thus not obtained. 

6 Summary 

In summary, the electronic and magnetic structure of doped 
and compressed NiS2 was discussed in great detail within 
different approximation in the framework of DFT. In par- 
ticular, LDA/GGA, GGA+C/, and a hybrid functional 
composed of 20% EECE mixed to LDA have been com- 
pared to Hartree-Fock and GW calculations. The metal- 
insulator transition with applying pressure or doping in 
NiS2 can be described properly only when using the hy- 
brid functional. The GGA+C/ fails to describe correlated 
metals since is works only well in case of integer orbital 
occupation and long range order. However, as the electron- 
lattice interaction is very strong in NiS2, and the results 
very sensitive to changes in the lattice structure, the re- 
sults of LDA/GGA calculations already allows to iden- 
tify the microscopic origin of the metal-insulator tran- 
sition, namely the dependence on crystal field splitting 
and bonding-antibonding splitting. Additionl corrections 
to LDA/GGA are then necessary to provide the insulating 
band structure of NiS2. 

Therefore, the magnetic properties of NiS2, on the 
other hand, are described best within GGA. LDA under- 
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estimates the magnetic order, thus no magnetic ground 
state is found. The employed hybrid functional overesti- 
mates the magnetic order, even when only a small amount 
of Fock exchange is mixed to LDA. The GGA+U approach 
cannot cpature the magnetic properties as expected for 
metallic compounds away from half band filling. 

Thus, the Fock exchange applied to the transition metal 
Ni in the EECE approach gives reasonable results as com- 
pared to experimental data as far as the electronic proper- 
ties are concerned as well as in the insulating as the metal- 
lic phases of NiS2 . The ordering of the magnetic moments 
is highly overestimated. The magnetic order present in the 
EECE approach alters the band structure in comparison 
to Hartree-Fock. 
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Electronic and magnetic properties of NiS2-a;Se^: a comparative study 

Cosima Schuster 
Institut fiir Physik, Universitdt Augsburg, 86135 Augsburg, Germany 
(Dated: November 17, 2009) 

The metal-insulator transition and the problem of d-electron delocahzation are investigated in 
the pyrite system NiS2-a;Sea; under pressure using density functional theory (DFT). We test several 
approximations, including the generalized gradient approximation (GGA), the GGA-|-i7 approach, 
and hybrid functionals. In addition we apply the GW approximation and perform Hartree-Fock 
calculations. The important role of the chalcogen dimers in the electronic structure is discussed 
within GGA, which sufficiently describes the role of the lattice in the metal-insulator transition. 
In addition, the magnetic phase diagram is determined. However, the electronic properties are 
inadequately described - the insulating ground state of NiS2 cannot be obtained - and the magnetic 
order is slightly overestimated. If correlations are taken into account within GGA-\-U , the insulator 
is found, but the non-magnetic ground state of the doped samples is not accessible. Mixing Fock 
exchange with local approximations as GGA correctly reproduce the insulating ground state in 
NiS2, and metallic ground states with doping and pressure. The insulator is - in contrast to earlier 
suggestions - of Mott type, i. e. the gap opens between Ni 3det, states. Due to the Fock term, the 
magnetic order is strongly overestimated. The applied GW calculations are not able to correct the 
metallic character of NiS2 and underestimate the gap by one order of magnitude. 

PACS numbers: 73.20.-r,74.25. Jb,85.25.Am 

Keywords: Metal-insulator transition, density functional theory, magnetic pliase diagram 



I. INTRODUCTION 

Crystal field splitting, electron-electron interaction, 
and the id band width are almost equal in the pyrites 
MX2 (M=Fe,Co,Ni,Cu,Zn; X=S,Se) which gives rise to 
a large variety of electrical, magnetic and optical prop- 
erties in these compounds. Of special interest is the 
metal-insulator transition in NiS2 with Se-doping Jj] or 
under pressure • The transition is not accompanied by 
a change in the lattice symmetry and is commonly be- 
lieved to be driven by the electron-electron interaction. 
In addition, a low-temperature antiferromagnetic phase 
is observed both in the insulating and metallic regime 
of NiS2-a;Sea; [ij and by applying pressure in NiS2 
The electronic and optical properties of NiS2 were first 
investigated already 30 years ago [J, and re-investigated 
experimentally in the late nineties [1, 0, 0, H, • The re- 
newed interest in transition-metal oxides and their metal- 
insulator transitions traces back to the attempt to under- 
stand the phase diagram of the high-Tc superconducting 
perovskites. Moreover, the attempt to use the metal- 
insulator transition of V2O3 or VO2 for technical appli- 
cations triggered many studies up to now [l0| . 

Theoretical attempts using density functional theory 
and the local density approximation (LDA) to describe 
the electronic structure of the pyrites date back to 1987 
pT| , with the focus on FeS2 , where the important role of 
the S(e)-S(e) dimers present in the pyrite structure has 
been pointed out. The role of electron-lattice interaction 
when applying pressure was studied ten years later, again 
in FeS2 |12l |. Other pyrites are mentioned in [ll|, and 
their electronic structure is explained by filling the band 
structure with the c?-electrons of the transition metal by 
going from Fe {M^^s^ ) to Zn (3di°4s2). The band insu- 
lators are described quite well within the local density ap- 



proximation, however the microscopic origin of the insu- 
lating nature of NiS2 is not well understood. Since more 
advanced approximations within DFT are available now, 
several attempts to improve the material specific theory 
have been performed. Recently, Perucchi et al. [13| have 
discussed two different mechanism for the metal-insulator 
transitions based on LDA calculations. Compression and 
expansion of the lattice alter the metallic state, albeit 
through two different microscopic mechanisms. On the 
other hand, Kunes et al. [ll|, using LDA-I-DMFT, have 
traced back the metal-insulator transition to a single con- 
trol parameter, namely the size of the p gap. They con- 
firmed, however, that compression and doping influences 
this control parameter in different ways. Both studies 
concentrate on the electronic structure in the paramag- 
netic phase and neglect the magnetism. 



Hence, we clarify the role of the S-S dimers and dis- 
cuss the energy scales mentioned in flp . [l^ . For the 
calculations we use the wien2k code '15| and the gener- 
alized gradient approximation (GGA) in the parameter- 
ization of PBE 16]. In a next step we determine the 
magnetic phase diagram. Sec. III. Surprisingly, GGA de- 
scribes the magnetic phase diagram quite well, slightly 
overestimating the phase boundaries as compared to ex- 
perimental data. In contrast to GGA, LDA always fa- 
vors the non-magnetic ground state. Calculations using 
improved exchange-correlations functionals are presented 
in Sec. IV. We present resuhs of GGA-f ;7(SIC) [H and 
EECE hybrid [H [H calculations. To compare with we 
have performed also GW calculations on top of LDA us- 
ing ABINIT [2O1, Jlj] . We focus on the advantages and 
shortcomings of the different techniques with regard to 
correlated materials. 
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II. THE ROLE OF THE DIMERS FOR THE 
ELECTRONIC STRUCTURE 

The pyrite structure is a simple cubic lattice with space 
group Pa3. It is best described in terms of the NaCl 
structure with the transition metal in one sub-lattice and 
the center of mass of the chalcogenide pairs in the other. 
Thus, the S/Se atoms form distorted octahedra around 
the Ni atoms. In the following, we use the experimen- 
tal lattice constant of NiS2, a = 5.687 A ^2^. Applying 
hydrostatic pressure is simulated by a shortening of the 
lattice constant - 5 %/10 % volume reduction correspond 
to 1.7 %/3.5 % reduction of the lattice constant and to 
5/12 GPa, respectively. In the attempt to understand 
the phase diagram with doping we concentrate on the 
end points with x = Q and x = 2. In addition we con- 
sider also a; = 1, i. e. NiSSe. For NiSe2 we also use the 
experimental lattice constant, a = 5.960 A [l^. In case 
of NiSSe, we study two different configurations within 
the pyrite unit cell - recent STM measurement on doped 
samples show no sign of a superstructure [l^l - and use 
the theoretically obtained volume, atheo = 5.87 A which 
is about 1 % larger than estimated in experiments [l|, [l^l . 
The first configuration contains S2 and Se2 dimers, the 
second only S-Se pairs. Thereby the second configura- 
tion is favored by 0.015 Ryd when the volume is fixed. In 
general, S-S, S-Se, and Se-Se pairs are present in doped 
samples Q. 

Using the mentioned lattice constants the internal po- 
sitions are obtained by a force minimization. The dis- 
tances are summarized in Tab. HI In NiS2, the S2 distance 
is found to be about 2.1 A, slightly (2 %) larger than 
the with experimental data d = 2.07 |Tl[. Applying 
moderate pressure, the dimer decreases only slightly from 
2.104 to 2.095, thus the dimer-distance is changed by 
0.4% when changing the lattice constant by 3.5%. Apply- 
ing even more pressure, corresponding to a 20% volume 
effect, it decreases only to 2.07 A. Also when expanding 
the volume drastically, the dimer-distance vary hardly, 
compare [ij] . For this reason, the metal- insulator transi- 
tion applying pressure is related to a volume effect, since 
the dimer-distance is nearly less constant when varying 
the volume. Turning to the doped samples, we find a 
large variation of the dimer-distances with doping. In 
NiSe2, the Se-Se distance is given hy d = 2.49 A, which 
is about 4% larger than the experimental 11]. Worth 
mentioning is that the Se-Se distance is larger than the 
Ni-Se distance. Turning to NiSSe, we observe that in 
NiSSe I the the S-S/Se-Se distances are comparable to 
dimer-distances of the pristine materials, see Tab. reftab. 
In NiSSe II, on the other hand, the S-Se distance corre- 
sponds to the averaged value of d—2.sA.. 

Next, we study the volume/dimer-distance effect on 
the electronic structure, see Fig. [U where we show the 
density of states (DOS) of NiS2. Generally, the presence 
of S/Se2 dimers in the pyrites leads to a sphtting of the 
bonding and anti-bonding S 3s states as well as S 3p 
states, as indicated. A change of the volume is expected 



System 


a (A) 


d: Ni-S/Se 


d: S(e)-S(e) 


AEs (eV) AEp (eV) 


NiS2 


5.28 


2.21 


2.07 


2.2 


11.0 


NiS2 


5.49 


2.31 


2.10 


3.2 


9.9 


NiS2 


5.69 (exp.) 


2.40 


2.10 


3.4 


9.0 


NiS2 


6.04 


2.56 


2.11 


2.3 


7.2 


M^CO^ T 

iNiobe i 




2.48 (a) 




3.35 (S) 


8.6 (S) 






2.44 (Se) 


2.42 


2.3 (Se) 


8.0 (Se) 


NiSSe II 


5.87 


2.41 (S) 


2.30 


1.9 


7.7 






2.50 (Se) 








NiSe2 


5.96 


2.48 


2.49 


1.8 


7.4 



TABLE 1: Lattice constants and typical bond lengths - given 
in A- of the systems under consideration, NiS2, compressed 
NiS2, expanded NiS2, NiSSe, and NiSe2. In addition, the 
energy differences of the bonding and anti-bonding splitting 
of the chalcogen (S, Se) s and p states are given. 
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FIG. 1: (Color online) Density of states of NiS2. We show the 
S 3s, S 3p3, and Ni 3d states. The bonding and antibonding 
parts of the S 3p states are indicated. 



to modify the bandwidth. With increasing volume the 
bandwidth should decrease, and vice versa. However, 
the situation in NiS2 is more complicated as pointed out 
for FeS2 in [l2|. Following the detailed analysis in [l2| 
the peaks at —16 and —12 eV are assigned to S 3s states. 
These states are observed in the XPS valence spectra at 
— 15 and —11 eV [6]. Thus, the splitting is correctly de- 
scribed in GGA, but the position shifted by 1 eV. Domi- 
nant S 3p contributions are found at the lower band edge 
(at —7 eV) and the upper band edge (at 2 eV). They 
form in itself a double peak structure, compare [I4I. The 
Ni 3d eg states - located between —4.5 eV and —2.54 eV 
and between — 1 eV to 1 eV - overlap with S 3p states. 
The Ni 3d t2g states, located at —2 eV, show almost no 
overlap with the S 3p-states. 

Applying pressure, see Fig. [21 upper panel, where we 
show a comparison of the DOS of compressed and un- 
compressed NiS2, a rigid shift of the t2g by —0.5 eV due 
to the enhanced crystal field splitting is observed in the 



DOS. Due to this shift the Ni t2g- and the S p-states 
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FIG. 2: (Color online) Upper panel: Density of states of NiS2 
at the experimental lattice constant and compressed NiS2. 
The straight (red) line corresponds to zero pressure, the dot- 
ted (green) line to 12 GPa. Lower panel: Density of states of 
NiS2 (red line) in comparison to the density of states of NiSe2 
(green line). 



overlap. The states below —3 eV are subject to an addi- 
tional shift, so that the lower band edge is about 0.6 eV 
lower for a = 5.49 Athan for a — 5.69 A. Hence the over- 
all band width is broader and the bonding anti-bonding 
splitting increases when compressing the lattice. In ad- 
dition to the shift, the binding S 3p-states show a slight 
broadening, but the rigid band shift is the dominant ef- 
fect. 

Upon doping, see Fig. [21 lower panel, where we show a 
a comparison of the DOS of NiS2 and NiSe2 , the bonding- 
antibonding splitting is reduced due to the increased Se- 
Se distance, both for the s and p states, see Tab. U as 
compared to NiS2. In NiSSe I, the low energy double 
peak structure due to the s states splits into four peaks 



due to the presence of both S-S and Se-Se dimers with 
different dimer lengths. This could be an experimental 
signature whether a structure with both S-S and Se-Se 
dimers is realized or only S-Se pairs are present. The 
dip at Ep is more pronounced in the NiSSe II than in 
the NiSSe I structure. In addition, the conduction band 
width is smaller in structure II (taken at the same vol- 
ume). To summarize, the dominant effect of doping on 
the DOS/band width is the decrease of the splitting of 
the S 3p-states and the resulting stronger overlap of the 
Ni 3d with the S p states, as commonly suggested. The 
stronger hybridization is identified as the driving mech- 
anism of the metal-insulator transition upon doping. 

Up to now, we have analyzed the density of states 
over the whole band width. Turning to the underlying 
band structure, we take closer look at the states near the 
Fermi level. Interestingly, the changes in the dispersion 
e(fc) around Ep are similar when applying pressure and 
with doping. For comparison, we plot the according band 
structures in Fig. [3l The changes in the band structure 
near the Fermi level with pressure can be focused on two 
symmetry points: The band at F (0.3 eV) moves down- 
wards, almost touching or crossing the Fermi level, re- 
spectively. Second, the four degenerate bands at M are 
shifted towards the Fermi level. In NiSe2, the band at R 
moves from the Fermi level to about 0.4 eV above. The 
hybridization due to band shifts or decreased bonding- 
antibonding splitting does not manifest itself near Ep. 



III. MAGNETIC PHASE DIAGRAM 

The phase boundaries of the AFM ordering of type I 
with pressure and doping can be reproduced within the 
DFT calculations using GGA. From the experience with 
LDA/GGA studies of other magnetic materials it is ex- 
pected that LDA underestimates and GGA overestimates 
magnetic order and magnetic moments, see for example 
(25| . In NiS2 , we consider an antiferromagnetic alignment 
of type A, i. e. the magnetic structure consists of ferro- 
magnetic planes which are coupled antiferromagnetically. 
Using this structure no supercell is needed, in accordance 
with Q. The magnetic properties - calculated within 
GGA - are summarized in Tab. [Ill The given magnetic 
structure is stable within the GGA calculations, frustra- 
tion effects play a minor role. As expected, NiS2 turns 
out to be a high-spin compound with /i = 0.73/iB. With 
increasing pressure, the magnetic moment as well as the 
energy difference to the non-magnetic case decreases. In- 
creasing the pressure up to 12 GPa, the antiferromagnetic 
ground state is still favored, but the energy difference to 
the non-magnetic solution tends to zero. NiSe2 turns out 
to be non-magnetic, with zero moment. At the critical 
doping (NiSSe), we find also an antiferromagnetic order, 
which can be suppressed by applying pressure. Unfortu- 
nately, a lower boundary of the magnetic transition by 
using LDA could not be determined since LDA favors 
the non-magnetic solution in all cases. This can be re- 
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FIG. 3: (Color online) Band structure near Ep of - from top 
to bottom - NiS2, NiS2 with 10% volume reduction and for 
NiSez. 



lated to a volume effect. Within LDA the ground state 
volumes are much smaller than the experimental or GGA 
volumes which may have an impact on the spin-polarized 
calculations. For completeness, the data for the energy 
differences using LDA are also given in Tab. HH 

The electronic properties change due to the magnetic 
order. Band structure and density of states are heavily 
influenced by the magnetism, see Fig. HI where we show 
DOS and band structure of antiferromagnetic NiS2. The 
degeneracy of the bands is lifted from a fourfold to a 
twofold degeneracy. States near Ep, especially at R and 
M, are emptied. As a result, the dip in the DOS at the 
Fermi level becomes more pronounced. 
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FIG. 4: (Color online) Band structure and density of states 
of NiS2 with anti-ferromagnetic spin configuration. 



compound 

NiSa 

NiSa 

NiS2 

NiSz 

NiSa 

NiSSe 

NiSSe 

NiSSe 

NiSea 



TABLE II: 

within LDA/GGA calculations for the systems under investi- 



Vxc a (in 


A) mag. order moment (in 


/is) AS in eV 


GGA 5.69 


afm 


0.73 


2.1 


LDA 5.69 


non 


0.20 


-0.14 


GGA 5.49 


non 


0.45 


0.17 


LDA 5.49 


non 


0.11 


-0.17 


GGA 5.28 


non 







GGA 5.69 


afm 


0.4 


0.034 


GGA 5.87 


afm 


0.5 - 0.6 


0.88 


LDA 5.87 


non 


0.1/0.2 


-0.14 


GGA 5.96 


non 







Magnetic ; 


ground state and 


magnetic 


moment 



gation. Ai5 = (£„ 



-pol " 



-i?mag)/(number of Ni) is the energy 



gain due to the magnetic ordering per Ni atom and is given 
in meV. 



IV. GGA+f/, FOCK EXCHANGE 

Local approximations to DFT like LDA and GGA fail 
to describe the insulator [13|. Recently, a variety of in- 
sulating transition-metal oxides have been investigated 
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using LDA+C/ LDA+DMFT Fock exchange 

or hybrid functional [l^, or the GW approximation 
[2^ . In this work, we apply GGA+C/, hybrid functionals, 
Hartree-Fock, and GW calculations to NiS2. We use the 
implementation of GGA+U, comparing the fully local- 
ized limit (SIC) as applied to NiO [13] and the around 
mean-field treatment (AFM) for the double count- 
ing correctiomand the EECE treatment fl9i| within the 
wien2k code [15| . 

To start with, we perform constraint GGA calculation 
to obtain a reasonable value for U. Following the treat- 
ment in j28| and assuming a Ni d^ configuration as in 
NiO, we obtain U^g = U - J = 6.39 eV for NiSj. f/eff 
decreases to 5.89 eV for NiSSe, and further to 4.94 eV 
for NiSe2. The different chemistry such leads to a con- 
siderable decrease of the correlation energy with doping. 
With decreasing volume, UeS is constant. In the follow- 
ing, we always use the antiferromagnetic spin structure as 
discussed in the previous section. The DOS, calculated 
with Uos = 6.4 eV and the SIC correction, is plotted 
in Fig. [5l upper panel, and shows a rather poor agree- 
ment with experimental data 0] ■ The Ni 3d states are 
shifted to too low energies, forming a triple peak struc- 
ture at the lower band edge. The S 3p states are homo- 
geneously spread over the whole energy range. Using the 
experimental value of U = 5 eV the gap values, namely 
A^^C = 0.6 eV and A^'^'^ = 0.3 eV, are in the range of 
the experimental values, but the Ni 3d states are located 
at —4 eV instead of —2 eV. A similar behavior was found 
in case of the ferromagnetic metallic pyrite C0S2 i29| . 

While GGA+U is successful for insulators, the descrip- 
tion of correlated metals where the requirement of an in- 
teger orbital occupation and well established long range 
order [30| is no longer fulfilled may fail. In case of NiS2 
under pressure, GGA+U correctly shows that this sys- 
tem is a non-magnetic metal. The results for NiSe2 on 
the other hand, are very sensitive to the choice of the 
method. Using SIC the ground state is a ferromagnetic 
metal, in case of AFM the antiferromagnetic metal is 
nearly degenerate to a ferromagnetic metal with two dif- 
ferent moments. 

In a next step we perform calculation using Fock ex- 
change for the Ni within the mufhn tin spheres (exact 
exchange for correlated electrons, EECE) and LDA for 
the correlation potential. This approach is known to give 
quite reasonable results for correlated electron systems 
[isl [19! . Instead of a parameter U and the choice for the 
double counting corrections which give deviating results 
for insulating NiS2 and metallic NiSe2, the mixing pa- 
rameter of the hybrid functional has to be determined. 
As expected 100 % EECE overestimates the gap by a fac- 
tor of 5. It also shifts the Ni 3d states to the lower band 
edge, i. e. to -10 to -8 eV. In this way, 100 % EECE is 
comparable to U = 10.7 eV. The gap is induced within 
the S p states, and not as expected between the Ni 3d 
(Mott-Hubbard type) or Ni 3d and S 3p states (charge 
transfer type). If 25 % EECE is added to LDA the Ni 
3dt2g states are found at —4 eV. The DOS shape corre- 




-8 -6 



-4 -2 

E-EF(eV) 



FIG. 5: (Color online) Total density of states of NiS2 (red 
line) and contributions of the Ni 3d (green and blue line) and 
S 3p (pink line) orbitals: we compare the GGA+U calculation 
(left-hand side) with the hybrid functional with 20% EECE 
(right-hand side). 



sponds to the DOS obtained for U — 5.4 eV (suggested 
U value from |11| ) . Decreasing the amount of EECE fur- 
ther, the gap closes below 20%; therefore we perform the 
following calculation with a hybrid functional where 20% 
EECE is mixed to LDA. Then, gap in the DOS opens 
within the Ni-S hybridized states, located between —2 
and 1 eV, at the preformed dip in the GGA calculation, 
as can be seen in Fig. [51 The bonding-antibonding S 3p 
states dominate at the lower valence and higher conduc- 
tion band edge. The double peak structure of the empty 
states observed in XAS can be traced back to a Ni 3deg 
double peak and the S 3p antibonding double peak. Us- 
ing this hybrid functional, an antiferromagnetic solution 
is obtained, although LDA itself favors the non-magnetic 
ground state. The calculated moments become smaller 
under pressure and with doping but do not vanish. 



6 



The gap of NiS2 is very sensitive to the S-S dimer- 
distance d. Using a structure with d = 2.15 Athe band 
gap amounts to 0.3 eV, whereas it is 0.6 eV with d = 2.10 
A. Increasing the dimer distance to d — 2.36 A, i.e. using 
the structural parameter of NiSe2 but the lattice constant 
of NiS2, the gap closes. A closer look at the weighted 
band structure shows that the gap is located at F; the 
lower band has dominant Ni character, whereas the upper 
band has dominant S character. The S character of this 
band shows up only at F, compare FeS2 [l2|- According 
to the band structure the metal-insulator transition un- 
der pressure is triggered by the band with Ni character 
which provides the states at the Fermi level between F 
and X. With doping the metallic state is attributed to the 
states between X and M. In this way, a metallic ground 
state is obtained for NiSe2 due to the increased dimer- 
distance. On the other hand, reducing the volume also 
closes the gap at constant dimer-distance. The results for 
the band structure are shown in Fig. [51 where we com- 
pare insulating NiS2 with metallic NiS2 under pressure 
and metallic NiSe2. 

To summarize, GGA+U is not able to describe the sys- 
tem properly, neither with the calculated nor estimated 
U values. Hybrid functionals using Fock exchange can 
reproduce the insulator and give a reasonable DOS. This 
approach fails, however, in describing the magnetic phase 
diagram by overestimating the ordered phase. Thus, we 
have identified a system where the advantage of hybrid 
fimctionals versus GGA+U can be clearly demonstrated. 



V. QUASI-PARTICLE PICTURE: HF AND GW 

Corrections to the energy eigenvalues in DFT can be 
estimated by first order many-body perturbation theory 
with respect to the difference in exchange correlation self- 
energy and exchange correlation potential of DFT. De- 
scribing the dielectric properties of semiconductors, the 
GW approximation has become the method of choice 
[3l| . Thereby, the non-local, energy-dependent self- 
energy is approximated by the Green's function and the 
screened Coulomb interaction, S = GW [s^l- The screen- 
ing is calculated on the level of the random phase approx- 
imation. Within the "one-shot" procedure Go and Wo 
are calculated from LDA/GGA, but this GqWq approach 
gives appropriate results only if GW is a small cor- 
rection to LDA. Transition-metal oxides can be treated 
rather well using self-consistent approaches. For exam- 
ple, the electronic structure of VO2 is described correctly 
[26j . Even properties of V2O3 can be described [11]. In 
NiS2, the electron-lattice coupling is strong and found 
to be the driving mechanism of the metal-insulator tran- 
sition, thus it can be expected that the GW approach 
is able to provide a description of the insulating state. 
Similar to VO2, in the case of NiS2, a self consistent 
scheme has to be applied to get rid of the metallic start- 
ing point. In the following we use the GW implemen- 
tation of the ABINIT package [2l|. We compare, see 



NiS2 0.2EECE 

3.0 -I , 




R FX M r 



NiSe2 

3.0 -| 




R r X M r 



FIG. 6: Band structure obtained using a hybrid functional 
with 20% EECE. Top: NiS2; middle: NiS2 under pressure; 
bottom: NiSe2. 

[3^ . Hartree-Fock calculations with self-consistent, but 
static GW schemes. The self-energy can be written as 
a sum of the exchange self-energy (static screened ex- 
change), the correlation self-energy (Coulomb- hole term) 
and the dynamical (energy-dependent) self-energy, and 
the approximations used within GW are classified along 
this partitioning. Here, we compare the screened ex- 
change with the screened exchange plus Coulomb-hole. 
Another static GW app roximation is the quasi-particle 
self-consistent GW [35|. As we start from the Kohn- 
Sham orbitals obtained within in LDA, spin-polarized 
calculations are in accordance with non-polarized calcu- 
lations since LDA favors a non-magnetic ground state. 
For simplicity we chose the k-points where the correc- 
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method spin 




k-point 


gap 


HF non- 


■pol 


r 


0.3 eV 


HF pol 




r 


0.4 eV 


scrEX pol 




M 


0.175 eV 


COHSEX non- 


pol 


M 


0.021 eV 


COHSEX pol 




M 


0.025 eV 


QPscGW pol 




M 


0.027 eV 



TABLE III: Energy gaps obtained within different approxima- 
tion used for the self consistency, namely screened exchange, 
screened exchange plus Coulomb hole, and model GM^[3^. 



tions are calculated out of the optimal k-set. The given 
k-points Kl (1/8,1/8,1/8) and K4 (3/8,3/8,3/8) are lo- 
cated between T and R, and K2 (3/8,1/8,1/8) and K3 
(3/8,3/8,1/8) between X and M. As it can be assumed 
that the correction do not vary in the vicinity of the cho- 
sen k-points, Kl should give the correction at F, K2 at 
X, K3 at M, and K4 at R. 

On all GW levels, gaps are induced to the band struc- 
ture of NiS2 at the chosen k-points, where the smallest 
gap occurs at K3 (M). The results are summarized in 
Tab.Hni 

The behavior is in accordance with the GGA calcula- 
tions as we can assign the dip at Ep of the GGA DOS to 
the bands crossing at M. However, it is found that one 
band crosses the Fermi level between R and F. For this 
reason, GW calculations cannot describe the insulating 
state of NiS'2, even though a gap is formed at all given 
k-points. Therefore, calculations of compressed NiS2 and 
NiSe2 lead to very similar results for the gaps. Moreover, 
the metallic phases show larger gaps at these k-points, 
but the crossing is present in each case. 

On the other hand, the smallest gap is loacted at F 
within Hartrec-Fock, in accordance with the results dis- 
cussed in the previous section. In addition, the gap values 
resulting from the calculations with 20 % EECE and from 
the Hartree-Fock calculation are astonishingly of similar 
magnitude. However, the total density of states shows no 
gap at Ep but a large gap between —5 to —1 eV, despite 
the direct gaps at the chosen fc-points. 

As a major shortcoming, the GW and Hartree-Fock 
calculations, on top of the non-magnetic starting point 
generated by the LDA, do not generate a magnetic 
ground state as the EECE does. The local magnetic min- 
imum is thus not obtained. 



VI. SUMMARY 

In summary, the electronic and magnetic structure of 
doped and compressed NiS2 was discussed in great de- 



tail within different approximation in the framework of 
DFT. In particular, LDA/GGA, GGA+U, and a hybrid 
functional composed of 20% EECE mixed to LDA have 
been compared to Hartree-Fock and GW calculations. 
The metal-insulator transition with applying pressure or 
doping in NiS2 can be described properly only when us- 
ing the hybrid functional. The GGA+U fails to describe 
correlated metals since is works only well in case of in- 
teger orbital occupation and long range order. However, 
as the electron-lattice interaction is very strong in NiS2, 
and the results very sensitive to changes in the lattice 
structure, the results of LDA/GGA calculations already 
allows to identify the microscopic origin of the metal- 
insulator transition, namely the dependence on crystal 
field splitting and bonding-antibonding splitting. Ad- 
ditionl corrections to LDA/GGA are then necessary to 
provide the insulating band structure of NiS2. 

Therefore, the magnetic properties of NiS2, on the 
other hand, are described best within GGA. LDA under- 
estimates the magnetic order, thus no magnetic ground 
state is found. The employed hybrid functional over- 
estimates the magnetic order, even when only a small 
amount of Fock exchange is mixed to LDA. The GGA+U 
approach cannot cpature the magnetic properties as ex- 
pected for metallic compounds away from half band fill- 
ing. 

Thus, the Fock exchange applied to the transition 
metal Ni in the EECE approach gives reasonable results 
as compared to experimental data as far as the electronic 
properties are concerned as well as in the insulating as 
the metallic phases of NiS2 . The ordering of the magnetic 
moments is highly overestimated. The magnetic order 
present in the EECE approach alters the band structure 
in comparison to Hartree-Fock. 
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